Gamma-ray constraints on the infrared 
background excess 



M. Mapelli a , R. Salvaterra b , A. Ferrara a 

l SISS 'A /International School for Advanced Studies, Via Beirut 4, 34014 Trieste, 

Italy 

h Universita delVInsubria, Via Valleggio 11, 22100 Corao, Italy 



Abstract 

Motivated by the idea that the recently detected near-infrared (1.2-4 |im) excess 
over the contribution of known galaxies is due to redshifted light from the first 
cosmic stars (Salvaterra & Ferrara 2003), we have used the effect caused by photon- 
photon absorption on gamma-ray spectra of blazars to put constraints on extra- 
galactic background light (EBL) from the optical to the far-IR bands. Our analysis 
is mainly based on the blazar H 1426+428, for which we assume a power-law unab- 
sorbed spectrum. We find that an EBL model with no excess over known galaxies 
in the near-infrared background (NIRB) is in agreement with all the considered 
blazars; however, it implies a very peculiar intrinsic spectrum for H 1426+428. Ad- 
ditional data on the blazars 1ES1101-232, H 2356-309 and PKS 2155-304 exclude 
the existence of a strong NIRB excess consistent with Kelsall's model of zodiacal 
light subtraction (ZL); the COBE/DIRBE measurements, after Wright's model ZL 
subtraction, represent a firm NIRB upper limit. The constraints on the optical EBL 
are weaker, due to the fact that predictions from different optical EBL models are 
often comparable to the experimental errors. In the mid-infrared the SPITZER mea- 
surement of vl v =2.1 nW m~ 2 sr -1 at 24 fim, gives a good fit for all the considered 
blazars. 
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1 Introduction 



The study of the Extragalactic Background Light (hereafter EBL) might pro- 
vide unique information to understand many crucial astrophysical questions, 
including, among others, the first cosmic star formation (Salvaterra & Ferrara 
2003), the evolution of galaxies (Totani & Takeuchi 2002; Totani et al. 2001), 
the role of dust emission (Granato et al. 2001). Yet, the EBL determination in 
the optical and infrared bands poses challenging observational and theoretical 
difficulties, often leading to discrepant interpretations. For example, Madau & 
Pozzetti (2000) suggest that the optical EBL can be completely explained by 
galactic counts, while Bernstein et al. (2002) claim the presence of an optical 
excess. Moreover, the presence of an excess (Matsumoto et al. 2000; Wright 
2001) in the near infrared background (hereafter NIRB) crucially depends on 
the subtraction of the zodiacal light (Kelsall et al. 1998; Wright 1997, 1998; 
Wright & Reese 2000). For the middle and far infrared background (hereafter 
MIRB and FIRB) various models have been proposed (Totani & Takeuchi 
2002; Silva et al. 2004; Primack et al. 1999) which cannot be discriminated on 
the basis of currently available data (Elbaz et al. 2002; Metcalfe et al. 2003, 
and references here). Finally, very strong constraints to the MIRB can be 
provided by the SPITZER 1 satellite (Papovich et al. 2004), which indicates 
a background total flux of 2.7^oi 11 W m 2 sr^l at 24 /xm. In the next years 
SPITZER (and other projects like ASTRO-F and SPICA) will likely be able 
to reconstruct the complete map of the MIRB and FIRB. 
Waiting for the upcoming infrared facilities, an alternative way to put inde- 
pendent constraints to the EBL is to resort to the so-called photon-photon 
absorption of the high energy tail of the blazar spectra. It is known that 
gamma rays in the GeV-TeV energy bands can be absorbed by softer (mainly 
optical and infrared) photons, via electron-positron pair production (Nikishov 
1962). Stecker, de Jager & Salamon (1992) pioneered the method on the spec- 
trum of 3C 279 to first study the EBL. Their attempt stimulated a plethora of 
works on the subject (Stecker, de Jager & Salamon 1993; Madau & Phinney 
1996; Stecker & de Jager 1997; Aharonian et al. 2003; Konopelko et al. 2003; 
Stecker 2003; Costamante et al. 2004; Dwek & Krennrich 2005) that were en- 
couraged by the availability of new measurements of blazar spectra in the TeV 
regime based on the Cherenkov imaging technique (Aharonian et al. 1999b; 
Krennrich et al. 1999 and Aharonian et al. 2002b for Mkn 421; Aharonian et 
al. 1999a, 2001a, 2002b for Mkn 501; Aharonian et al. 2003 for H 1426+428). 

In this paper we calculate the optical depth due to photon-photon absorption 
using the most recent optical and infrared background measurements (Bern- 
stein et al. 2002; Matsumoto et al. 2000; Papovich et al. 2004), combined with 
state-of-art theoretical models of the NIRB (Salvaterra & Ferrara 2003) and of 
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the MIRB-FIRB (Totani & Takeuchi 2002). We compare these results with the 
Cherenkov imaging data available for six blazars (H 1426+428, Mkn 421, Mkn 
501, PKS 2155-304, 1ES1101-232, H 2356-309), to derive some constraints on 
the EBL. The final aim of the study is to assess if the most recent EBL data 
and models exclude or require the presence of an excess with respect to galaxy 
counts in the optical and in the near-infrared regime. In addition, we try to 
test the predictions of the proposed EBL models, especially in the optical, 
near and middle infrared ranges. 

The paper is organized as follows. Sec. 2 presents the method adopted to de- 
rive the optical depth for the photon-photon absorption and to compare it 
with blazar data. In Sec. 3 we briefly review the most important EBL data 
and the EBL models considered. Sec. 4 is devoted to the description of the 
Cherenkov data available for the three blazars that we analyze. In Sec. 5 we 
discuss the results and Sec. 6 summarizes our conclusions. 



2 The photon-photon absorption 



The optical depth due to photon-photon absorption can be written (Stecker 
et al. 1992; Madau & Phinney 1996) as 

r{E)= Jdz—Jdx [ ~ X) J den(e)<r(e,E,x), (1) 

-1 Hh 



where dl/dz is the proper line element 2 



dl 



[(l + z)£(z)}-\ 



dz H 
c is the speed of light and 



£{ Z ) = n M (i + zf + n A + (i - n A - n M )(i + zf 



1/2 



In eq. 1, x = cos#, 8 being the angle between the directions of the two 
interacting photons. E = E (1 + z) is the observed energy of the blazar 
photon and e = e (l + z) is the observed energy of the background photon; 
z em is the redshift of the considered blazar; finally, n(e) is the specific number 



We adopt the following cosmological parameters: Hubble constant Hq=71 km 
s- 1 Mpc- 1 , ft M =0.27, ft A =0.73, which are in agreement with the recent WMAP 
determination (Spergel et al. 2003). 
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density of background photons. The energy threshold for the interaction, e t h, 
is defined by: 

_ 2m^ 

eth ~ E{l-x) (2) 



where m e is the electron mass. 

The photon-photon absorption cross section (Heitler 1960; Gould & Schreder 
1967) is given by 



2/3(/3 2 -2) + (3- / 3 4 )ln 



.1-0, 



(3) 



where cft is the Thomson cross section and 



2m?c 4 

1 - 



Ee(l-x) 



1/2 



The cross section a peaks sharply at A ~ hcE/(2m 2 e c 4 ) ~ 2.4(E/TeV)/xm 
(where A is the wavelength of the interacting background photon). Fig. 1 
shows the behavior of cr(A) integrated over the angle 9. Fig. 2 shows how the 
cross section depends on the interaction angle 9. 

Substituting into eq. 1 the relation n(e) = n(e )(l + z) 3 (Madau & Phinney 
1996), we finally obtain: 



2 em xo 1 

C f [ \ -\- z) f f 
T ^ = 2~W J dZ £(z) ) dx ( l ~ x ) J de n(e )cr(e ,Eo,x,z)] (4) 

^ -1 e th0 

n(e ) is related to the observable quantity F(e ), the background photon flux 
at redshift z = 0, by the simple relation: 

n(e ) = 47T 7T- cm" 3 erg" 1 (5) 



We have calculated numerically t(E) using a three-dimensional integral based 
on the method of Gaussian quadratures (Press et al. 1992). The integration 
accuracy varies between 3 x 10" 5 % and 8 x 10" 3 % in the range 0.1-10 TeV. 
t(E) has been estimated for different values of F(e ), whose choice has been 
guided by both observational data and theoretical models (see next Section). 



4 




Fig. 1. Photon-photon absorption cross section, cr(A), as a function of the back- 
ground photon wavelength, A, integrated over the interaction angle, for three values 
of the observed energy of the blazar photon: 100 GeV (dashed line), 1 TeV (dotted) 
and 10 TeV (solid). &t = 6.652 xlO -25 cm 2 is the Thomson cross section. 




Fig. 2. Isocontours of the cross section (in units of o~t) as a function of 9 and A for 
E=l TeV. 



5 



To discriminate among these different estimates of t(E) we need a compar- 
ison with observational data of blazar spectra, which we have performed in 
two independent ways: 

(1) The first method assumes a theoretical model of the unabsorbed blazar 
spectrum. In particular, we adopt a simple power-law spectrum 
(dN/dE) unabsorbed oc E~ a for the blazar H 1426+428, the main target of 
our analysis. For the other blazars considered (Mkn 421 and Mkn 501) 
the unabsorbed spectrum is better fitted by a power-law with an expo- 
nential cut-off (Konopelko et al. 2003). Some theoretical models suggest 
that the blazar spectrum must present a cut-off in the TeV range (Inoue 
and Takahara 1996; Tavecchio, Maraschi & Ghisellini 1998; Fossati et al. 
2000). 

The absorbed spectrum is then obtained by convolving the unabsorbed 
spectrum with t(E), 



and by changing the spectral index a to obtain the best fit to the observed 
blazar spectrum. 

(2) The alternate method consists in inverting equation 6. In this case we 
apply t(E) directly to the observational data, to derive the unabsorbed 
spectrum and check if the intrinsic spectral index a is consistent with the 
current theoretical predictions. This procedure does not require any a 
priori assumption about the unabsorbed spectrum shape. An important 
sanity check is to verify whether the values of a derived from the two 
methods are consistent. 



3 The optical-infrared background 

3. 1 Observations of the EBL 

The major difficulty in measuring the EBL arises from the subtraction of the 
interplanetary dust scattered sunlight (i.e. zodiacal light, ZL) contribution. 
This problem presents different aspects and complexity according to the range 
of wavelengths observed. A vast amount of literature is present on the subject 
(see Hauser & Dwek 2001 for a complete review); here, we only concentrate 
on the the EBL measurements and models used in our work. 




(6) 
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3.1.1 Galaxy counts 



Deep optical and near infrared galaxy counts give an estimate of the EBL frac- 
tion coming from normal galaxies. Madau & Pozzetti (2000) derived the contri- 
bution of known galaxies in the U BVIJHK bands from the Southern Hubble 
Deep Field imaging survey. In particular for the U, V, B and I bands (corre- 
sponding to the wavelengths A = 3600, 4500, 6700 and 8100 A) they found a 
mean flux respectively 2.87lg;l, 4.57tg;$, 6.74tJ;jjl and 8.04l£;| in units of 
10~ 6 erg s _1 cm~ 2 sr _1 . 



3.1.2 Optical excess 

Estimates of the optical EBL based on photometric scans across dark nebulae 
(Mattila 1976; Spinrad & Stone 1979) and on photoelectric measurements 
(Dube et al. 1977, 1979) provide upper limits that are higher than the flux 
given by galaxy counts alone at the same wavelengths. In particular the most 
recent work by Dube et al. (1979) provides an upper limit of 2.6 x 10~ 5 erg 
s _1 cm~ 2 sr _1 at A= 5115 A. Bernstein et al. (2002) measured the mean flux 
of the optical EBL at 3000, 5500, and 8000 A, using the Wide Field Planetary 
Camera 2 (WFPC2) and the Faint Object Spectrograph, both on board the 
Hubble Space Telescope, combined with the du Pont 2.5 m Telescope at the Las 
Campanas Observatory. They found for these three band a mean flux of the 
EBL respectively 12.0±^ 3 7 , 14.9t$;ij, and 17.6± 2 |jj in units of 10" 6 erg s" 1 cm" 2 
sr -1 , considerably higher than the contribution of the galaxy counts alone. In 
this band not only ZL is likely to provide a substantial contribution, but 
also terrestrial airglow, and dust-scattered Galactic starlight (diffuse Galactic 
light) might represent a potential problem for the measurement. The impact 
of such systematic errors has led Mattila (2003) to question the claim by 
Bernstein et al. of the discovery of an optical EBL excess. 



3.1.3 Near Infrared Background: DIRBE and NIRS data 

The available NIRB data come from the Diffuse Infrared Background Exper- 
iment (DIRBE) on board of the Cosmic Background Explorer (COBE) and 
from the Near InfraRed Spectrometer (NIRS) on board of the InfraRed Tele- 
scope in Space (IRTS). The DIRBE instrument provided a survey of the sky in 
10 photometric bands at 1.25, 2.2, 3.5, 4.9, 12, 25, 60, 100, 140, and 240 fim us- 
ing a 0.7° x 0.7° field of view. A summary of the DIRBE results can be found 
in Hauser et al. (1998). The NIRS instrument covers the wavelength range 
from 1.4 to 4.0 /zm with a spectral resolution of 0.13 /xm. Matsumoto et al. 
(2000) made a preliminary analysis of the NIRS data, estimating the NIRB on 
the basis of the 5 NIRS observation days unperturbed by atmospheric, lunar 
and nuclear radiation effects. 
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Both the DIRBE and the NIRS data show an excess in the NIRB with respect 
to galaxy counts. An estimate of this excess depends on a critical point: the 
subtraction of the contribution of the ZL from the measurements. There are at 
least two models of ZL. The model described in Kelsall et al. (1998) exploits 
the temporal variability of the signal caused by looking at the sky through 
different amounts of the interplanetary dust as the Earth orbits the Sun. 
Wright & Reese (2000) noted that the high value of the EBL flux derived by 
Kelsall from the DIRBE data at 25 jum might indicate that a residual ZL flux 
could remain after subtraction. Thus they suggest a different approach, which 
requires that the EBL signal at 25 /im after ZL subtraction is zero (Wright 
1997, 1998). In practice, subtracting the ZL both with the Kelsall and with 
the Wright method the presence of a NIRB excess is unquestionable, even if 
the amount of this excess depends on the ZL model assumed. 



3.1.4 Mid Infrared Background: SPITZER data 



Before the SPITZER satellite, the only MIRB available data came from the 
ISOCAM deep extragalactic surveys. From the analysis of the ISOCAM num- 
ber counts Elbaz et al. (2002) computed an EBL flux (integrated down to 50 
/xJy) of 2.4 ±0.5 nW m~ 2 sr -1 at 15yum (68% confidence level). In deriving this 
value, Elbaz et al. took into account, among other surveys, of 15 /im counts 
from a portion of the ISO gravitational lensing survey. Metcalfe et al. (2003) 
use a full lensing survey (covering Abell 2218, Abell 2390 and Abell 370), and, 
by integrating from 30 /zJy up-wards, obtain an EBL flux of 2.7 ± 0.62 nW 
m -2 gr -i at i5^ m (68% confidence level). This result is consistent with the 
upper limit of 5 nW m~ 2 sr _1 on the 15 /zm EBL, estimated by Stanev & 
Franceschini (1998). This upper limit is calculated from photon-photon ab- 
sorption effects on the spectra of Mkn 501; note, though, that this analysis is 
based on relatively poor-quality, old data (Aharonian et al. 1997). 



Recently, integrating to 60 /xJy the counts from the Multiband Imaging Pho- 
tometer on board of the SPITZER satellite (MIPS, Rieke et al. 2004), Pa- 
povich et al. (2004) found a lower limit to the EBL flux at 24/im of 1.9 ± 0.6 
nW m~ 2 sr -1 . Extrapolating to fainter flux densities, they derive an estimate 
of the total 24 fxm background of 2.71^7 n W m~ 2 sr -1 , in good agreement 
with the result of Metcalfe et al. (2003). Other lower limits to the EBL flux 
at wavelengths ranging from 3 to 10 /im, derived by SPITZER measurements 
(Fazio et al. 2004), are reported in Fig. 3. 
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3.2 Theoretical models of the EBL 



The somewhat sparse experimental information available on the extended 
wavelength range of the EBL forces us to resort to theoretical modeling in 
order to fully reconstruct its spectrum. In the following we summarize the 
main features of the EBL models used in this work. 

3.2.1 Contribution of Pop III stars to the NIRB 

The most plausible explanation of the NIRB excess is that the EBL in this 
wavelength range is due to the redshifted UV and optical light emitted by 
Pop III stars (Bond et al. 1986; Santos et al. 2002; Salvaterra & Ferrara 2003). 
In particular Salvaterra & Ferrara (2003) developed a model of the NIRB 
which, accounting for the most recent predictions of Pop III stellar spectra 
(Schaerer 2002) and IMF, nebular emission (i.e. the radiation coming from the 
nebula surrounding the star), and L ya photons scattered by the intergalactic 
medium, is able to fit the NIRS data (Matsumoto et al. 2000) and the DIRBE 
data with both the methods of ZL subtraction. Their best fit predicts, for a 
star formation efficiency /* = 0.1 — 0.5, depending on the adopted IMF, a 
transition from (very massive) Pop III to Pop II stars occurring at z ~ 9. This 
model is supported also by the analysis of the Infrared Background fluctuations 
performed by Magliocchetti et al. (2003). Such interpretation of the NIRB in 
terms of Pop III stars, although very intriguing, might be somewhat extreme 
in terms of the high star formation efficiencies and production of intermediate 
mass black holes, as pointed out also by Madau & Silk (2005). 
To test this model and to understand the physical nature of the transition, a 
better knowledge of the EBL at A < 1.2 /im is crucial (i.e. a confirmation or 
rejection of the optical excess measured by Bernstein et al.). 

3.2.2 Mid- and Far Infrared Background 

For the middle and the far infrared background a plethora of different models 
is present in the literature. Among the most complete ones, that by Totani 
& Takeuchi (2002) is also consistent with the ISOCAM and the SPITZER 
measurements. These authors construct a model of the near, middle and far 
infrared background based on the backward approach for the luminosity evolu- 
tion of galaxies (i.e. they infer the star formation history from the present-day 
galaxy optical-IR SEDs and chemical properties). Their model is characterized 
by a realistic treatment of dust and a physical determination of its tempera- 
ture. 

For our work we have considered also other models of the IR background, 
such as those presented by Primack et al. (1999) or by Silva et al. (2004). The 
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Table 1 

Summary of the considered EBL models. 



Model Optical Background NIRB 



CI Madau & Pozzetti (2000) 

MK1 Madau & Pozzetti (2000) Matsumoto et al. (2000) (K) a 

MK2 Bernstein et al. (2002) Matsumoto et al. (2000) (K) 

DW1 Madau & Pozzetti (2000) Wright (2001) (W) a 

DW2 Bernstein et al. (2002) Wright (2001) (W) 
a (K) and (W) indicate the ZL subtraction obtained using Kelsall's model and 
Wright's model, respectively. For the MIRB and FIRB we adopted always the 
Totani & Takeuchi (2002) model rescaled to the Spitzer data ( Papovich et al. 
2004). 

results obtained using these models, however, do not give as good fits to both 
the Spitzer and TeV data as that of Totani & Takeuchi. 



3.3 Summary of the adopted EBL model 



The results of our study (presented in Sec. 5) are obtained adopting the fol- 
lowing assumptions on the EBL. 

• For the optical background (i.e. between 0.3 and 1.2 /im) we have consid- 
ered both the values obtained from galaxy counts only, following Madau & 
Pozzetti (2000), and the case of a background excess (Bernstein et al. 2002 
data, including their upper and lower limits). 

• For the NIRB (1.2 < A < 4. /im) we have used both the NIRS data (Mat- 
sumoto et al. 2000; Salvaterra & Ferrara 2003) with the Kelsall's model of 
ZL subtraction, and the DIRBE data with the Wright model of ZL sub- 
traction (Gorjian et al. 2000; Wright & Reese 2000; Wright 2001; Wright k 
Johnson 2001). 

• For the MIRB and FIRB (A > 4. jim) we adopt the Totani & Takeuchi 
(2002) model, rescaled by a factor 1.2 in the 8-30 /im range to match the 
SPITZER 24 /im data point; in this wavelength range the uncertain contri- 
bution of spiral galaxies allows such rescaling. 

Finally we have considered a case in which the optical/NIRB in 0.3 < A < 
4 /im is contributed purely by galaxies as given by their counts (i.e. no excess) 
leaving the MIRB/FIRB as above. A summary of the considered models is 
given in Table 1. Fig. 3 shows these models and the considered data. In our 
work we neglect a possible redshift evolution of the EBL. This assumption 
is reasonable given the low redshift of the considered blazars, including the 
most distant one (z = 0.129). We also neglect the possible self-absorption 
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0.1 1 10 100 1000 

A (/xm) 



Fig. 3. EBL data and corresponding models. In both the two panels are shown: 
data from Bernstein et al. (filled squares), Madau & Pozzetti (filled circles), Elbaz 
et al. (cross), Metcalfe et al. (open triangle), Fazio et al. (open squares from 3 to 
10 fim) and Papovich et al. (open square at 24 fim). The dotted line indicates CI 
model. Upper panel: DW1 (dashed line) and DW2 (solid line) models. The DW 
data (Gorjian et al. 2000; Wright & Reese 2000; Wright 2001; Wright & Johnson 
2001) are represented with open circles. For comparison, the dot-dashed line shows 
the P1.0 model of Aharonian et al. 2005b. Lower panel: MK1 (dashed line) and 
MK2 (solid line) models. The MK data (Matsumoto et al. 2000) are represented 
with filled triangles. 
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of the blazar, i.e. the contribution to the photon-photon absorption given by 
infrared photons produced by the blazar, as it can be shown to be irrelevant 
(Protheroe & Biermann 1997). 



4 TeV blazar spectra 

Sufficiently high resolution blazar spectra in the TeV regime, obtained with 
imaging Cherenkov techniques, are available for at least three blazars: H 1426+428, 
Mkn 421 and Mkn 501. Data have become recently available for other three 
blazars: PKS 2155-304 (Aharonian et al. 2005a), 1ES1 101-232 and H 2356-309 
(Aharonian et al. 2005b). H 1426+428 is at relatively high redshift (z = 0.129); 
its TeV spectrum determination is therefore significantly less accurate than 
that for Mkn 421 and Mkn 501. For this blazar available data from CAT 
(Djannati-Atai et al. 2002), Whipple (Horan et al. 2002) and HEGRA (Aha- 
ronian et al. 2003) exist. These data must be distinguished in two sets. The 
first set includes the HEGRA data taken in 1999-2000, CAT data taken in 
1998-2000 and Whipple data taken in 2001. The second set is represented by 
HEGRA data taken in 2002 and it is characterized by a much lower flux level 
than for the previous campaigns. Mkn 421 and Mkn 501 are maybe the best 
observed blazars in the high energy gamma-ray band and are characterized by 
nearly equal redshifts [z = 0.031 and z = 0.034, respectively), considerably 
lower than that of H 1426+428. Because of this, differences in their spectra 
(for example in the cut-off energies) cannot be explained by different amounts 
of intervening absorption, but as due to intrinsic spectral characteristics (Aha- 
ronian et al. 2002b; Konopelko et al. 2003). The spectrum of Mkn 501 was 
measured up to 22 TeV both by HEGRA (Aharonian et al. 1999a, 2001a) and 
by Whipple (Samuelson et al. 1998; Krennrich et al. 1999): the observations 
of these two Cherenkov telescopes are in good agreement. Mkn 421 was also 
measured up to 20 TeV both by HEGRA (Aharonian et al. 1999b, 2002b) and 
Whipple (Krennrich et al. 1999, 2002), yielding similar fluxes. HESS spectra 
have been collected for PKS 2155-304 (z = 0.117), 1ES1101-232 (z = 0.186) 
and H 2356-309 (z = 0.165) (Aharonian et al. 2005a, 2005b). In Sec. 5 we will 
present spectral fits obtained for these blazars through the photon-photon 
absorption calculation, dwelling mainly on H 1426+428. 



5 Results 

As discussed in section 3.3, we have considered several models of optical and 
near infrared background (see Tab. 1). For all these models we have calculated 
the optical depth due to photon-photon absorption using eq. 4, and applied 
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Table 2 

X 2 values and best fit spectral indexes (a) of the considered models. 



CI 


MK1 


MK2 


DW1 


DW2 


X 2 8.66 


17.88 


19.86 


6.43 


7.57 


a 2.65 


1.60 


1.70 


1.80 


1.90 



Statistical analysis based on 12 observational data of the spectrum of H 1426+428 
(reported in Fig. 7) with 2 parameters (the spectral index a and a normalization 
factor) . 




0.1 1 10 

E 7 (TeV) 

Fig. 4. The solid line represents the best fit to the spectrum of H 1426+428 obtained 
for the model CI (see Table 1). The shaded area indicates the ±1<7 uncertainty in- 
troduced by the error in galaxy counts on the absorbed spectrum. The observational 
data reported here are from CAT 1998-2000 (filled squares), Whipple 2001 (filled 
triangles), HEGRA 1999-2000 (filled circles), HEGRA 2002 (open circles). All error 
bars are at 1 a (Aharonian et al. 2003). 



it to the theoretical unabsorbed blazar spectrum (eq. 6), changing the spec- 
tral index and the normalization until the best fit was found through the x 2 
method. Let us consider now the results obtained for each different model. 
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5.1 Galaxy counts only 



We start from the analysis of the most conservative case (CI) in which the 
optical and the NIRB come only from normal galaxies as derived by Madau 
k Pozzetti (2000); while the MIRB/FIRB are from the rescaled Totani & 
Takeuchi model. Fig. 4 shows the best fit, obtained through x 2 minimization, 
to the H 1426+428 data for model CI. The best fit is obtained for spectral 
index a = 2.65 with x 2 ~ 9 (Table 2), derived considering 12 observational 
data (CAT 1998-2000, Whipple 2001, HEGRA 1999-2000) and 2 free param- 
eters (spectral index a and normalization factor). This spectral index seems 
to be uncommonly high for this blazar. In particular, for H 1426+428 previ- 
ous literature indicates a spectral index a = 1.9 (Aharonian et al. 2002a) or 
a — 1.5 (Aharonian et al. 2003). If we consider a spectral index a < 2, we 
obtain x 2 > 21, a considerably higher value. 

Another problem of the CI model is that the HEGRA point at E~ 5TeV is 
more than 2a away from the curve. Moreover, the best fit drops at E> 3 TeV, 
hence providing a very poor match to the observed shape of the absorbed 
spectrum derived from experiments, which flattens in the energy range 1-6 
TeV. In conclusion, although model CI cannot be rejected on the basis of x 2 
analysis alone, we consider it unlikely given its poor performance in terms of 
spectral slope and shape. 



5.2 Including a NIRB excess 

We next consider models including a NIRB excess in the range 1.2 — 4 /im 
as the Matsumoto et al. (2000) data with the Kelsall's subtraction of the 
ZL (hereafter MK), and the DIRBE data with the Wright subtraction of the 
ZL (DW); again, we fix the MIRB/FIRB according to the rescaled Totani & 
Takeuchi model in both cases. For the optical background we experimented 
with all available measurements (Madau & Pozzetti 2000; Bernstein et al. 
2002; Mattila 2003). 

Fig. 5 shows the photon-photon absorption optical depth for a blazar at 
redshift 0.129 (i.e. the redshift of H 1426+428), assuming the EBL models 
MK1+MK2 (left panel) and DW1+DW2 (right panel). The optical depth in 
the MK1 and MK2 cases (those which are characterized by the MK data in 
the NIRB) is considerably higher, especially at 1 TeV (where the absorption 
comes primarily from NIRB photons, see Fig. 1). Fig. 6 shows the best fits for 
the H 1426+428 spectrum, again with the NIRB from MK1 and MK2 models 
(left panel) and from DW1 and DW2 models (right panel). DW1 and DW2 
provide an excellent fit, with x 2 respectively 6.4 and 7.6 (Table 2). None of 
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Fig. 5. Optical depth for photon-photon absorption with z em = 0.129, the redshift 
of H 1426+428. The EBL models assumed are the following. Left panel: MK1 
{dashed line) and MK2 (solid); Right panel: DW1 (dashed line) and DW2 (solid 
line). Shaded areas refer to ±1<7 errors. 

the observational data is more than la away from the fit, that correctly re- 
produces the shape of the absorbed spectrum suggested by the data, with the 
plateau at E> 1 TeV. 

On the contrary, the MK1 and the MK2 cases are only marginally consistent 
with at least 4 Cherenkov data (CAT at ~ 1 TeV and ~ 0.6 TeV, HEGRA at 
~ 0.75 TeV and Whipple at ~ 0.78 TeV). Also the \ 2 is considerably higher 
than in the case of DW1 and DW2, even if not so high that we can reject 
MK1 and MK2 models. The good result obtained adopting DW1 and DW2 
models supports the evidence for a NIRB excess, although not as pronounced 
as suggested by the MK1 and MK2 models. 



5.3 Including the optical excess 

We finally consider the differences among various measurements of the optical 
background. The left and the right panel of Fig. 6 show the differences in the 
spectrum of H 1426+428 due to the assumption of an optical background due 
to galaxies only or including an "optical excess" as measured by Bernstein et 
al. In both cases the fit seems to favor an optical background consistent with 
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galaxy counts only (dashed line), although the differences between various 
types of optical background are smaller than the error bars in the blazar data. 
Thus more solid conclusions on the optical excess have to await for more 
precise data. 



5.4 The spectrum of H 1426+428 



From the two best models including the NIRB excess shown in Fig. 6, we can 
derive the shape of the intrinsic (unabsorbed) H 1426+428 spectrum. In both 
cases we find that a simple power-law of the form dN/dE cx E~ a represents 
well the unabsorbed spectrum. In particular we find that the best fit is given 
by a = 1.6 for the MK1 model, a = 1.7 for MK2, a = 1.8 for DW1 and a = 1.9 
for DW2. These values are consistent with previous works for H 1426+428, 
which estimate a slope for the unabsorbed spectrum of this blazar of a = 1.9 
(Aharonian et al. 2002a) or a = 1.5 (Aharonian et al. 2003). 
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Fig. 7. Fit of the H 1426+428 spectrum. We have assumed an EBL model given 
by: (a) optical background (0.3-1.2 fim): Madau & Pozzetti 2000; (b) near infrared 
background (1.2-4. /an): DIRBE (Wright 2001) with the ZL model of Wright & 
Reese (2000); (c) middle and far infrared background : model of Totani & Takeuchi 
(2002). The model of Totani & Takeuchi for the MIRB and the FIRB has been 
rescaled in the 8-30 fim range assuming different values of the EBL at 24 /xm, 
in particular assuming: 2.7 (solid line), 3.0 (dotted line), 3.3 (short dashed line), 
3.5 (dot-dashed line) and 3.8 (long dashed line) nW mT 2 sr _1 . The observational 
data reported here are from CAT 1998-2000 (filled squares), Whipple 2001 (filled 
triangles), HEGRA 1999-2000 (filled circles). 

5.5 Constraints on the SPITZER measurement 



So far we have always used the EBL value measured by SPITZER of 2.7 
nW m~ 2 sr _1 at 24 /im. It is interesting to check the consequences of an 
higher value, corresponding to the measurement upper limit of 3.8 nW m -2 
sr _1 . In Fig. 7 we assume an optical background given only by galaxy counts 
and a DIRBE NIRB (with the ZL model performed by Wright), which can 
be considered the fiducial model according to the previous analysis. For the 
MIRB/FIRB we adopt the Totani & Takeuchi model rescaled at 24 /zm to 
flux values in the range 2.7-3.8 nW m~ 2 sr _1 . Fig. 7 illustrates the effect 
of an increase of the 24 //m flux: it becomes more and more difficult to fit 
the Cherenkov data, even with substantial changes of the spectral index a. 
On the contrary, the best fit to the blazar spectrum is undoubtedly obtained 
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Table 3 

Spectral indexes (a) derived using the least square fits for our models CI, DW1, 
MK1 and for model P1.0 of Aharonian et al. 2005b. 





H 1426+428 


PKS 2155-304 


1ES1101-232 


H 2356-309 


CI 


2.5 


2.7 


1.4 


2.0 


DW1 


1.8 


2.0 


0.4 


1.0 


MK1 


1.4 


1.6 


-0.1 


0.1 


P1.0 






-0.1 


0.7 



adopting the SPITZER best value (2.7 nW m 2 sr l ). This result is clearly 
model dependent (the Totani & Takeuchi model has been assumed). 

5.6 An alternative analysis 

As discussed in Sec. 2, there is a different way to tackle the data analysis. So 
far we have assumed a theoretical shape (i.e. a power-law) for the unabsorbed 
blazar spectrum and convolved it with the optical depth for photon-photon 
absorption. Alternatively, one can apply t{E) directly to the Cherenkov data 
by inverting eq. 6 and derive the unabsorbed blazar spectrum without making 
any a priori assumption about its shape. This second approach is particularly 
indicated when the blazar cannot be fitted by a simple power-law (for exam- 
ple Mkn 421 and Mkn 501). As a sanity check, we now use this alternative 
approach to show that the two methods yield a coherent picture. 

Fig. 8 shows the results of such attempt for H 1426+428. We find that the 
CI model yields an intrinsic spectrum which is inconsistent not only with 
a power-law, but also with an exponential cut-off, as a significant rise in the 
spectrum above 3 TeV is seen. Models MK1, and especially DW1, do not show 
this peculiar spectral rise at high energies. Unfortunately, though, because of 
the large experimental errors in H 1426+428 data, we cannot assess clearly 
whether the intrinsic spectrum of this blazar is better fit by a simple/broken 
power-law or a power-law with an exponential cut-off. If we assume a power- 
law intrinsic spectrum, we can derive the best fit spectral index with the 
weighted least square method. With this method we derive a = 1.8 for the 
case DW1, a = 1.4 for MK1 and a = 2.5 for CI, consistent with what we 
found minimizing the x 2 (Table 3). 

Hence, the results of this analysis, and in particular the anomalous rise in 
the spectrum above 3 TeV, tend to disfavor the CI model; the DW1 model is 
found again to give the best fit to the data. 

Fig. 9 and Fig. 10 show the same procedure applied to Mkn 421 and Mkn 501. 
The calculated optical depths for these two blazars in the range 0.7-2 TeV are 
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Fig. 8. Unabsorbed spectrum of H 1426+428 obtained by convolving the pho- 
ton-photon optical depth with the data in Fig. 7 and assuming EBL models CI 
(open squares), DW1 (open circles) and MK1 (open triangles). The lines are the 
least square fits for the unabsorbed spectrum power-law index: CI (dot-dashed line, 
a = 2.5), DW1 (dotted, a = 1.8) and MK1 (dashed, a = 1.4). The error bars are 
obtained from the la error bars of the blazar data. 



in agreement with the results by Konopelko et al. (see Fig. 2 of their paper); 
above 2 TeV our results are slightly different, due to the higher MIRB flux of 
the EBL model (Malkan & Stecker 2001; Stecker 2003) they assume. The data 
for these two blazars have significantly higher quality than those available for 
H 1426+428; unfortunately, because of their low redshift, the effect of photon- 
photon absorption on the spectrum of these two blazars is too weak to provide 
additional constraints on the EBL. 

Figs. 11-13 show the results of this procedure applied to the new data available 
for PKS 2155-304, 1ES1 101-232 and H 2356-309, respectively. In the case of 
PKS 2155-304 both the spectrum derived from the DW1 model and that 
derived from the MK1 model show a peculiar peak at energy ~ 1 TeV. This 
peak is mainly due to the fact that these EBL models are characterized by a 
high flux at 1-4 /im and a rapid decrease at wavelengths > 4 //m. This finding 
is at odd with the conclusions from the analysis of H 1426+428. This is even 
more puzzling if we note that the two blazars are nearly at the same redshift. 
Instead, the unabsorbed spectrum derived from the model CI is quite smooth. 
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Dwek, Krennrich & Arendt (2005) use the theoretical emission models of Chi- 
appetti et al. (1999), based on the EGRET data (Vestrand, Stacy & Sreekumar 
1995), to put a constraint on the hardness of the spectrum of PKS 2155-304. 
We have compared the results of our models with EGRET data and Chiap- 
petti et al. models (Fig. 11). In agreement with Dwek et al. (2005), we find 
that the model CI is perfectly consistent with Chiappetti et al. (1999) model. 
It is outside the scope of the present paper to discuss the details of synchrotron 
self-Compton models. Nevertheless, taking into account the errors on the least 
square fits (> 0.8 for both DW1 and MK1 models), on the EGRET measure- 
ments and the uncertainties of the Chiappetti et al. model, we find that the 
DW1 model is close to be acceptable, whereas the MK1 model has, probably, 
to be rejected. On the other hand, we also have to take into account that 
EGRET data and observations in the TeV range are separated by a consider- 
able time interval (about 10 years); this can be crucial, given the variability 
of blazar spectra (Ghisellini, personal communication). Thus, the comparison 
between EGRET and TeV data, although interesting in principle, must be 
taken with some care. 

In Fig. 12 the results for the blazar 1ES1101-232 are presented. Aharonian et 
al. (2005b) have shown that the EBL models with a high excess in the NIRB 
must be rejected, since they imply an unabsorbed spectrum so hard that 
it would be difficult to explain it within the standard hadronic or leptonic 
scenario (Aharonian 2001b). This conclusion is supported by our study: our 
model MK1 gives a spectral index a = —0.1 (Table 3), which is too hard to be 
allowed by current theoretical models. However, for the model DW1 we obtain 
a spectral index a = 0.4 (Table 3), and for model DW2, in which the optical 
EBL matches the data by Bernstein et al. (2002), we get a = 0.5. Taking 
into account that the error on the spectral index from the least square fits is 
large ( > 0.3 for all the considered models), we conclude that the unabsorbed 
spectra derived both from DW1 and DW2 models are hard, but consistent 
with the lower limit a = 0.6 predicted for 1ES1101-232 by synchrotron self- 
Compton models (Ghisellini, personal communication). Then, EBL models 
with Wright-subtracted-ZL NIRB excess cannot be rejected on the basis of the 
HESS data of 1ES1101-232. In Fig. 12 we show, for comparison, the model 
P1.0 of Aharonian et al. (2005b). This model matches the COBE/DIRBE 
measurements with the Wright model of ZL subtraction, as the DW1 model 
does, but it has a higher flux from 4 to 10 /xm (Fig. 3) with respect to the latter. 
It is important to note that the model P1.0 gives a spectral index a = —0.1, 
which is as hard as predicted by MK1. This points out the importance of the 
spectral region 4-10 /zm to understand the EBL. 

Finally, Fig. 13 shows the blazar H 2356-309. Also in this case, the spectral 
index for the MK1 model is very hard (a = 0.1, Table 3), whereas both CI 
and DW1 models give unabsorbed spectra of acceptable hardness. 
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Fig. 9. Unabsorbed spectrum of Mkn 421 obtained by applying the photon-photon 
absorption to the 2000/2001 HEGRA data (Aharonian et al. 2002b). The adopted 
EBL models are as in Fig. 8. 



6 Conclusions 



In this paper, we have explored the effects produced by photon-photon absorp- 
tion on blazar spectra to put constraints on extragalactic background light 
from the optical to the far infrared bands. We mostly focused on the high 
redshift blazar H 1426+428, for which we adopted a power-law unabsorbed 
spectrum. This might be a simplifying assumption, as some theoretical models 
indicate a more complex spectral shape (Inoue and Takahara 1996; Tavecchio, 
Maraschi & Ghisellini 1998; Fossati et al. 2000). On the other hand, the data 
for H 1426+428 are not sufficiently accurate to suggest the existence of a 
cut-off or to exclude a power-law spectrum. 

We have considered three different models of the NIRB, based respectively on 
galaxy counts (CI) and on the presence of an excess (DW and MK). In the 
case of H 1426+428 the DW model provides the best fit both using the method 
presented in Sec. 5.2, and the alternative one of Sec. 5.6. The MK model has 
a higher \ 2 but cannot be rejected. For the same reasons, CI model cannot be 
excluded either; however, in this case a very peculiar blazar intrinsic spectrum 
must be assumed (Fig. 8). We conclude that the presence of a NIRB excess 
with respect to galaxy counts, at the level given by the DW model, seems to 
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Fig. 10. Unabsorbed spectrum of Mkn 501 obtained by applying the photon-photon 
absorption to the 1997 HEGRA data (Aharonian et al. 2002b; filled squares). 

be required to fit the blazar spectrum. 

The very recent HESS data of the blazars PKS 2155-304 (Aharonian et al. 
2005a), 1ES1101-232 and H 2356-309 (Aharonian et al. 2005b) allow us to 
put stronger constraints to the NIRB. In particular, from the analysis of all 
these blazar, we conclude that the MK models, based on the Kelsall method 
of ZL subtraction, must be rejected, since they imply too hard blazar spectra. 
The spectra of PKS 2155-304 and 1ES1101-232 are marginally consistent with 
DW models, based on Wright ZL subtraction, which can be considered as an 
upper limit to the NIRB. The model CI, without NIRB excess, is favored by 
all these new data. These findings are quite different from those we derived for 
H 1426+428. Since the available data for H 1426+428 are old, new measure- 
ment of this blazar are eagerly required, to shed some light on this puzzling 
inconsistency. 

The derived constraints on the optical EBL are weaker, due to the fact that 
deviations between different optical EBL models are comparable to the exper- 
imental errors. The fit tends to favor models without an optical excess over 
galactic light, contrary to the result obtained by Bernstein et al. (2002). A 
more solid conclusion on the amplitude of the optical EBL has to await for 
the next generation of Cherenkov Telescopes as MAGIC, VERITAS, HESS, 
GLAST, or infrared satellites (GIBER). 
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Fig. 11. Unabsorbed spectrum of PKS 2155-304 obtained by applying the pho- 
ton-photon absorption to the HESS data (Aharonian et al. 2005a). The adopted 
EBL models are as in Fig. 8. The lines are the least square fits for the unabsorbed 
spectrum power-law index: CI (dot-dashed line, a = 2.7), DW1 (dotted, a = 2.0) 
and MK1 (dashed, a = 1.6). The error bars are obtained from the la error bars of 
the blazar data. The crosses are the EGRET data (Vestrand, Stacy & Sreekumar 
1995) and the thin short dashed line represents the best fit synchrotron self-Compton 
model of Chiappetti et al. (1999). 



Finally, in the mid-infrared the SPITZER measurement of ul u =2.7 nW m~ 2 
sr _1 at 24 /im, combined with the EBL model by Totani & Takeuchi (2002), 
allows us to obtain a good fit for all the blazars available. Again, a tremendous 
advance on the determination of the MIRB/FIRB is expected from SPITZER 
and the next generation of infrared satellites as ASTRO-F, and SPICA. 



In summary, recent measurements of blazars in the TeV range seem to ex- 
clude the existence of a strong NIRB excess consistent with Kelsall's model of 
ZL subtraction. The COBE/DIRBE measurements, after Wright's model ZL 
subtraction, represent a firm NIRB upper limit. 
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Fig. 12. Unabsorbed spectrum of 1ES1101-232 obtained by applying the pho- 
ton-photon absorption to the HESS data (Aharonian et al. 2005b). The adopted 
EBL models are: CI (open squares), DW1 (open circles) and MK1 (open triangles). 
In addition we show also the results obtained for the P1.0 model of Aharonian 
et al. 2005b (crosses). The lines are the least square fits for the unabsorbed spec- 
trum power-law index: CI (dot-dashed line, a = 1.4), DW1 (dotted, a = 0.4), MK1 
(dashed, a = —0.1) and P1.0 (solid, a = —0.1). The error bars are obtained from 
the la error bars of the blazar data. 
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